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^ ■ Abstract 

; . Using kicked differential equations of motion with derivatives of noninteger orders, we obtain gen- 

l> 
^D \ eralizations of the dissipative standard map. The main property of these generahzed maps, which are 

called fractional maps, is long-term memory. The memory effect in the fractional maps means that 
their present state of evolution depends on all past states with special forms of weights. Already a 
small deviation of the order of derivative from the integer value corresponding to the regular dissipa- 
tive standard map (small memory effects) leads to the qualitatively new behavior of the corresponding 
attractors. The fractional dissipative standard maps are used to demonstrate a new type of fractional 
attractors in the wide range of the fractional orders of derivatives. 



Discrete maps are widely used to study the general properties of dynamical 
systems. In those cases when they can be derived from differential equations, 
their analysis gives the exact properties of the corresponding systems. In this 
article we derive discrete maps (fractional maps) from the fractional differential 
equations, which correspond to the fractional generalizations of the dissipative 
standard map [1]. We demonstrate how the attractors of the fractional maps 
are different from the attractors of the dissipative standard map. 

I. INTRODUCTION 

There is a number of distinct areas of physics where basic problems can be reduced to the 
study of simple discrete maps. Discrete maps as substitutes of differential equations have 
been used to study evolution problems in Refs. 2|-[6|. They lead to a simpler formalism, 
which is particularly useful in simulations. The dissipative standard map is one of the most 
widely studied maps. In this paper we consider fractional generalizations of the dissipative 
standard map which are described by fractional differential equations [THSl]. 

The treatment of nonlinear dynamics in terms of discrete maps is a very important step in 
understanding the qualitative behavior of systems described by differential equations. The 
derivatives of noninteger orders are a natural generalization of the ordinary differentiation 
of the integer order. The fractional differentiation with respect to time is characterized 
by long-term memory effects which correspond to intrinsic dissipative processes in physical 
systems. The application of memory effects to discrete maps means that their present state 
of evolution depends on all past states |l0l-ll5l bj 



Discrete maps with memory can be derived (see Ref. 17|]) from equations of motion with 



fractional derivatives. In Ref. 



ITJ I a fractional generalization of the standard map has been 



derived from a fractional differential equation. A fractional generalization of the dissipative 



standard map was also suggested in 17|. Unfortunately, in that generalization a dissipation 
was introduced by the change of the variable pn -^ —bpn- The map equations were not 
directly connected with a fractional equation of motion. In this paper we propose two 
generalizations of the dissipative standard map. The first one is derived from a differential 
equation with fractional damped kicks. The second generalization of the dissipative standard 
map is derived from a fractional differential equation (kicks are not fractional) . A nonlinear 



system with fractional derivatives perturbed by a periodic force exhibits a new type of 
chaotic motion which can be called the fractional chaotic attractor 16|. Fractional discrete 
maps [l7| are used to study new types of attractors of fractional dynamics described by 
kicked fractional equations. In this paper some fractional differential equations of motion 
of kicked systems with friction are considered. Corresponding discrete maps with memory 
are derived from these equations. The fractional generalizations of the dissipative standard 
map are suggested and these maps are used in computer simulations. 



II. DISCRETE MAPS WITHOUT MEMORY 



In this section, a brief review of discrete maps is considered to fix notations and provide 



convenient references. For details, see Refs. |2|-|6| 



A. Standard map 



Let us consider the equation of motion 

oo 

X + K sin(x) 2^ ^ ( ^ 



n =0, 



(1) 



n=0 



in which perturbation is a periodic sequence of delta-function-type pulses (kicks) following 
with period T = 1, i^ is the amplitude of the pulses. This equation can be presented in the 
Hamiltonian form 

oo 

X = p, p + K sm{x)y^6it — n) = 0. (2) 

n=0 

..t is well-known that these equations can be represented (see, for example. Chapter 5 in Ref. 
Sj) in the form of discrete map 

Pn+1 =Pn- K sin(x„). (4) 

Equations ([3]) and (jl]) are called the standard map. This map is also called the Chirikov 
map [^. 



B. Dissipative standard map 

The dissipative standard map [l|, |l8|, ll9|] is 



Xn+l — Xn + /^5^n+l + ^, (5) 

F„+i = e-^[r„ + £sin(X„)], (6) 

where /x = (e'^ — l)/q- The dissipative standard map is also called the Zaslavsky map. Note 

that a shift Vt does not play an important role and it can be put to zero [Vt = 0). The 
dissipative standard map with f2 = can be represented by the equations 

Xn+l = Xn + Pn+l, (7) 

Pn+1 = -hPn - Z sin(X„). (8) 

For the parameters 

Z = -e/ie"'?, Pn = lJ.Yn, b = -e-^ (9) 

equations ([7]) and (|8]) give Eqs. (|5]) and (Q with fi = 0. 

For b = —1 and Z = K, we get the standard map which is described by Eqs. (|3]) and (j4]) 
with T = 1. 

Note that for large g — )■ oo (for small 6^0) Eqs. i^^ and (^ with Z = —K shrink to 



the one- dimensional sine map proposed by Arnold 20|, 



Xn+i = Xn + Ksm{Xn). (10) 

C. Kicked damped rotator map 

The equation of motion for a kicked damped rotator is 

oo 

x + qx = KG{x)^5{t-nT\ (11) 

n=0 

It is well known [S] that Eq. ( ITTl) gives the two-dimensional map 

1 - e"''^ 

Xn+l=Xn^ [pn + KG{Xn)\, (12) 

Pn+i = e-'i^[Pn + KG{xn)]. (13) 



This map is known as the kicked damped rotator map. The phase volume shrinks each time 
step by a factor exp(— g). The map is defined by two important parameters, dissipation 
constant q and force amphtude K. These equations can be rewritten in the form 

It is easy to see that these equations give the dissipative standard map (jTj) and ([H]) with 
fi = if we use 

X„ = Xn, Yn = Pn, s = K, T = 1, G{x) = sin(a;). 

This allows us to derive dissipative standard map ([5]) and ([6]) from the differential equation, 

oo 

X + qX = 6sm{X)J2^(t-n)- (14) 

71=0 

These equations give the discrete map defined by Eqs. ([5]) and ([6]) with fi = 0. 

III. FRACTIONAL STANDARD MAP AND DISSIPATION 

A fractional generalization of the differential equation ([1]) has been suggested in Ref. 



ITJ I . The discrete map which corresponds to the fractional equation of order 1 < a < 2 was 
derived. This map can be considered as a generalization of the standard map for the case 
1 < a < 2. 

We consider a fractional generalization of Eq. ([1]) in the form 

oo 

oD^x + Ksm{x)^6(t-n) =0, (1 < a < 2), (15) 

n=0 

where qD" is the Riemann-Liouville fractional derivative |7H9|, which is defined by 

Here we use the notation D^ = cP/dt"^, and o-^" is a fractional integration [7H9|. 
Defining the momentum as 

Pit) = oDr'xit), 



and using the initial conditions 

ioDr'x)iO+) = p,, ioDr'x)iO+) = b, (17) 

it is possible to derive the equation for the fractional standard map. 

Proposition 1. The fractional differential equation of the kicked system / I73]) is equiva- 
lent to the discrete map 

1 _"^ h 

Xn+i = —-Y,Pk+iVa{n -k + l) + ^ _ {n + ir-', (18) 

^ ^ k=0 ^ ' 

Pn+1 =Pn- Ksm{Xn), (1 < « < 2), (19) 

where the function Va{z) is defined by 

V^{z) = z"-' -{z- l)"-\ (20) 



Proof of this Proposition is given in Ref. 21 1 



A fractional generalization of the dissipative standard map suggested in Refs. |l|, [l8| can 

be defined by 

1 " 
x„+i = -— -Vpfc+iV„(n-A; + l), (l<a<2), (21) 

i (a) ^ — ^ 

^ ^ fc=0 

Pn+i = -bpn - Z sin(x„), (22) 

where the parameters are defined by conditions ([9]). For 6 = — 1 and Z = K Eqs. (!2T|) 
and fl22]) give the fractional standard map with T = 1. Note that this fractional dissipative 
standard map is not derived from a fractional differential equation. This map is derived 
by Pn ~^ ~bpn in the fractional standard map. Fractional dissipative standard map can be 
derived from fractional differential equations. In this paper, we derive two fractional gen- 
eralizations of the dissipative standard map which are obtained from fractional differential 
equations. 

IV. FRACTIONAL DERIVATIVE IN THE KICKED TERM AND THE FIRST 
FRACTIONAL DISSIPATIVE STANDARD MAP 

In this section we suggest the first fractional generalization of differential equation (TTTj) 
for a kicked damped rotator. In this generalization we introduce a fractional derivative in 

6 



the kicked damped term, i.e. the term of a periodic sequence of delta-function-type pulses 
(kicks), and derive the corresponding discrete map. 

Consider the fractional generalization of equation f llip in the form 



oo 

D',X{t) - qDlX{t) = esin (^Dfx) ^5(t -n), (0 < /3 < 1), 



(23) 



n=0 



n 



where g e M, and ^D^ is the Caputo fractional derivative [£] of the order < /3 < 1 defined 
by 

?DfX = ,ll-^DlX = —1- f -^^, (0 < /3 < 1). (24) 



r(l-/3)7o it-rf dr 
Here we use the notation D] = d/dt, and qI" is a fractional integration |7H9|. For /3 = 

fractional equation (l23l) gives equation (fTTl) . Note that we use the minus on the left-hand 

side of Eq. ( !23l) . where q can be a positive or negative value. Fractional derivative qD^X 

is presented in the kicked damped term. 



Proposition 2. The fractional differential equation of the kicked system [2B\] is equiva- 
lent to the discrete map 



X. 



n+l ^^n 



Xn + 



1 -e~'i 



^n+1) 



Yr 



n+l = e" 



n, 



Yn + e sin L ^ n+i 1^2-/3 (g, A; 

^ '^'^ fc=0 

where the functions W2-ii{a,h) are defined by 

W2-p{a, b) = a^-ie'*(^+^) [r(l - /3, ab) - 1(1 - (3, a{b + 1)) 
and r(a, 6) is the incomplete Gamma function 

/oo 
y^-'e-ydy. 



Proof. Fractional equation fl23|) can be presented in the Hamiltonian form 

X = Y, 

oo 

Y-qY = e sin( ^DfX) J] 5(t - n), 



n=0 



where < /3 < 0, and q eR. 
Between any two kicks. 



Y -qY = 0. 



(25) 
(26) 



(27) 



(28) 



(29) 



For t e {tn + 0, tn+i - 0), the solution of Eq. ([29]) is 



r(t„+i - 0) = r(t„ + 0)6". 

Let us use the notations t„ = nT, with T = 1 and 

X„ = X(tn - 0) = hmXfn - e), 

e-i-O 

V; = l^(tn-0) = hmF(n-e). 
For t G (tn — e, t„+i — e), the general solution of Eq. (128!) is 

oo „t 

r (t) = F„e''(*-*") + £ ^ sin( ?A1^) / dre'^^'-^'^Sir 

m=0 >^tn-e 

Then 

>;.+i=e''[r„ + £sin(o^Df„^) 
Using Eq. (1331) . the integration of the first equation of (l28l) gives 



(30) 



m 



Xn+l — Xn — 



F„ + 5sin(o^Df„^) 



(31) 



(32) 



(33) 



(34) 



Let us consider the Caputo fractional derivative from Eqs. (l33l) and (!M|) . It is defined by 
the equation 

1 /•*" c/r dX{T) 



Using l^('r) = dX{T)/dT, this relation can be rewritten as 



(0</3<l). 



\J Tn "p 



1 ^-1 /-i 



r-tfc+1 y(7-)(ir 



-t1/3^ 



(tn - r) 



(35) 



where t^+i = t^ + 1 = (fc + 1), and 4 = fc, such that to = 0. For r G (tfcjtfc+i), equations 
(150]) and dH]) give 



y(r) = r (tfe + 0)e«^"-*^^ = Y{tk+i - 0)e-V 



Qr.lij-tk) 



Then 






Y, 



fc+i 



tji ct 



gg(tn-ife+i-^)^-/3^^ ^ rfc+ie«(*"-*'=+i) 



tn— tfc + l 



tn—tk 



z-^e-'^'dz 



tn—tk + \ 



As a result, equation ( l36l) gives 



9(<n-tfc) 



(36) 



**+! Y{T)dT 

t, {in - ry 



Ffc+ig^-ie''^"-^-!) r(l - /3, g(t„ - t.+i)) - r(l - /3, g(t„ - 4)) 



(37) 



Here r(a, 6) is the incomplete Gamma function (1271) . where a and 6 are complex numbers. 
Using (l35l) and (1371) . we obtain 

n-l 

?^2^^ = f7r3mE^^'+i^2-/3(g,fc-n), (0</3<l), (38) 



fc=0 



where 



W^2-/3(a, 6) = aP-^e"^^^^^ r(l - /3, a6) - r(l - /3, a(6 + 1)) 
Substitution of Eq. (jHH]) into Eqs. ([33]) and (IHl) gives 



n-l 



K 



+i = e^ 



^n + gsin(^ ^ ^ Ffc+iW^2-/3(g, A; 



-n\ 



k=0 



Xn+1 — Xn 






n-l 



^n + £ sin(^— — — J2 Yk+iW2-^iq, k - 



n 



fc=0 



Equations ( HOj) and fHT]) can be presented in the form of Eqs. ( 125|) and (126|) . 
This ends the proof. D 



(39) 



(40) 



(41) 



The iteration equations f l25p and (12B]1 define a fractional generalization of the dissipative 
standard map. For /3 = this map gives the Zaslavsky map (E]) and (^ with 



/x=(l-e"'^)/g 



(42) 



and 17 = 0. 



V. FRACTIONAL DERIVATIVE IN THE UNKICKED TERMS AND THE SEC- 
OND FRACTIONAL DISSIPATIVE STANDARD MAP 

In this section we suggest a fractional generalization of the differential equation for a 
kicked damped rotator with fractional derivatives in the unkicked terms and derive the 
corresponding discrete map. 

We consider the fractional generalization of equation (TTT]) in the form 

oo 

oD'^Xit) - qoD^X{t) = £sin(X) 5^5(t - n), (43) 

n=0 

where 

qeR, l<a<2, 13 = a-l, 

and o-D" is the Riemann-Liouville fractional derivative |7H9|, which is defined by Eq. flT6l) . 
This equation has fractional derivatives in the unkicked terms, i.e. on the left-hand side of 
Eq. fH51) . We use the minus in the left-hand side of Eq. fH5|) . where q can have a positive 
or negative value. 

Proposition 3. The fractional differential equation of the kicked system [J^ is equiva- 
lent to the discrete map 

1 " 

^n+l = :p7 TT J^ ^+1 W^{q, k - U - 1) , (44) 



r„+i = e^[y„ + £sin(X„)J, (45) 

where the functions Wa{a,b) are defined by 

Wa{a, b) = a^~"e''(*+^) [r(a - 1, ab) - T{a - 1, a{b + l))j , (46) 

andr{a,b) is the incomplete Gamma function [27\ ). 

Proof. Let us define an auxiliary variable ^(t) such that 

?Dr"e=x(t), (47) 

where qD^~°' is the Caputo fractional derivative (^^. Using 

olt'' ^Dt"-^ = e(t) - m, (0 < 2 - a < 1), (48) 

10 



we obtain 



oD^X = Di ,It-^X = Di oin oDf--^ = Dfim - m) = A'^, 



(49) 



and 

oD^X = D] oll-^X = D] olt^'X = 

Substitution of Eqs. (HH]), (EU]), and (gT]) into Eq. ^ gives 



(50) 



DU - QDk = e sin( o^A'"°0 $^ ^(t - n), (1< a < 2). 



ra=0 



The fractional equation fl5T|l can be presented in the Hamiltonian form 



(51) 



^ = Y, 



Y -qY = esm{^Df-''^)J2Kt-n), (1< a < 2, qe 



n=0 



Using Eq. fl26l) of Proposition 2, we obtain 



n-l 



Y 



n+i = e^ 



Yn + e sin ( — — ^^ Yk+iWa{q, k - n) 

Vi a — 1 ^ — ^ 



For (X„, Yn), we use equation ( 138|) in the form 



fc=0 



n-l 



^ ^ k=0 



n) 



As a result, we have 



1 

Xn+i = T^T -T V n+iH^a(g, A: - n - 1), 

Fa — 1 -^^ 

^ ^ fc=0 



Y 



n+l 



F„ + e sin(X„ 



where Wa{a,b) is defined in Eq. ( H6|) . This ends the proof. D 



If we use the variables 



Pn = /iFn, 6 = -e", Z = -fxee'^ 



(52) 



(53) 
(54) 
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FIG. 1: Structures of the chaotic attractors for different values of K obtained after 10^ iterations 
(r = 5, f] = 0); a = 1.9975 in (a)-(e). (a) K = 12.83; (b) K = 12.87; (c) K = 12.93; (d) zoom of 
(c); (e) K = 12.97; (f) q = 2, K = 12.93. 



then equations (jH]) and fH5|) give 



X. 



n+l 



/^" 



— Y,Pu+iWM,k-n-l). 



(55) 



r(a-l) 

^ ^ fc=0 

P„+i = -6P„-Zsin(X0, (56) 

These equations can be considered as a fractional generahzation of the dissipative standard 
map equations ([7]) and ([8]) with f] = 0. For a = 2, this fractional dissipative standard map 
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FIG. 2: Structures of the fractional chaotic attractors for K = 12.93 and different values of a 
obtained after 10^ iterations (F = 5, il = 0). 

gives the dissipative standard map that is described by Eqs. ([7]) and ([8]). 

VI. NUMERICAL SIMULATIONS 

Numerical simulations were performed for the second fractional dissipative standard map 
(Eqs. fl55|) and f l56|) ). First we used our code to reproduce the results presented in Fig. 1 
from pj;] for the structures of the chaotic attractors of the dissipative standard map at the 



13 



window of the ballistic motion near K ^ An {q = —5, K = eexp(g), and used in Ref. 'l9| F 
is equal to —q) for the fractional standard map with a = 2 and obtained a perfect agreement 
(an example is given in Fig. 1(f)). As a decreases slightly from a = 2 to a = 1.9975, the 
window of the ballistic motion shrinks and moves to the higher values of K. Already for 
a = 1.9975 in Figs. l(a)-l(e) the ballistic motion appears for K > 12.86 and disappears 
at K = 12.97. The window is completely closed at a ~ 1.9969. The structures of two 
symmetric attractors with disjoint basins which appear within the window (Figs. 1(b) and 
1(c)) is also very different from the structures of the dying attractors of the dissipative 



standard map 



19| . The attractor in Fig. 1(d) evolves from period 8 trajectory to period 
4, period 2, and, finally, period 1 trajectory slowly moving in the direction of the upper left 
corner with the step of the order of 10^^. 

When a decreases further, the structures of the fractional chaotic attractors evolve in 
the manner presented in Fig. 2, where one can find one-scroll, two-scroll, and four-scroll 
fractional chaotic attractors, strongly deviating from the chaotic attractor of the dissipative 
standard map Fig. 1(f) (see also 19|). The problem of existence of multiscroll fractional 



chaotic attractors was considered in Ref. 221] but for the fractional differential equations 
with the Caputo derivatives. For values of a near 1 fractional chaotic attractor turns into 
period 2 and for smaller values period 1 attracting trajectories. Fig. 3. 

VII. CONCLUSION 

The suggested discrete maps with memory are generalizations of the dissipative standard 
map. These maps are potentially useful for description of fractional dynamics of complex 
physical systems. We consider viscoelastic and dielectric materials as examples of complex 
media (physical systems) whose dynamics could be described by the fractional dissipative 
standard maps and corresponding fractional differential equations. This assumption is based 
on the following. 

(1) It is well known that viscoelastic materials can be described by fractional differential 



equations (see, for example, Ref. 23|]). The fractional dissipative standard map could 
be employed to model the one-dimensional simplification of the equations of viscoelastic 
materials in which a perturbation is a periodic sequence of delta-function-type pulses (kicks) 
following with some period. 
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FIG. 3: Attracting trajectories for K = 12.93 and small values of a obtained after 10^ iterations 
(r = 5, J1 = 0). 

(2) For a wide class of dielectric materials the dielectric susceptibility follows a fractional 



24 



2511 over ex- 



power-law frequency dependence, which is called the "universal" response 
tended frequency ranges. The electromagnetic fields in such dielectric media are described 
by differential equations with fractional time derivatives 26|, \2j\. These fractional equa- 
tions for electromagnetic waves in dielectric media are common to a wide class of materials, 
regardless of the type of physical structure, or of the nature of the polarizing species. We 
assume that the fractional maps could be applied to dielectric media in which a perturbation 



15 



is a periodic sequence of kicks. 

The suggested fractional dissipative standard maps demonstrate a chaotic behavior with 
a new type of attractors. The interesting property of these fractional maps is long-term 
memory. As a result, a present state of evolution depends on all past states with the weight 
functions. The fractional dissipative standard maps are equivalent to the correspondent 
fractional kicked differential equations. Note that to derive discrete maps an approximation 
for fractional derivatives of these equations is not used. 

Computer simulations of the suggested discrete maps with memory prove that the non- 
linear dynamical systems, which are described by the equations with fractional derivatives, 
exhibit a new type of chaotic motion. This type of motion demonstrates a fractional gener- 
alization of attractors. 

that in the case q = (fractional standard map) the frac- 



It has been shown in Ref. 



tional discrete map demonstrates a new type of attractors such as slow converging and slow 
diverging trajectories, ballistic trajectories, fractal-like structures, and chaotic trajectories. 



At least one type of fractal-like sticky attractors in the chaotic sea can be observed 28 1 



for the fractional standard map. The properties of stability and existence of the fractional 

The 



attractors in the asymptotic sense for q = have also been described in Ref. p 
attractors presented in Figs. 1-3 are quite different from the corresponding regular chaotic 
attractors and g = attractors. The detail classification of these attractors and corre- 
sponding chaotic motion will be considered in the nearest future and published in the next 
paper. 
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